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Abstract 

The anisotropie eleetron transport in the (001) plane of sublimation-grown CI 2 -NDI (naphthalene di- 
imide) single erystals is analysed over a temperature range between 175 K and 300 K. Upon eooling from 
room temperature to 175 K the eleetron mobility along the direetion of preferred transport monotonously in- 
ereasesfrom 1.5 em^/Vsto2.8 em^A^s aeeordingto a distinet temperature relation of At first glanee, 

these eharaeteristies allude to a eoherent, i.e. band-like eharge earrier transport predominantly governed by 
inelastie seattering with aeeoustie phonons. However, as we will demonstrate, the experimental mobility 
data ean be eonsistently deseribed within the framework of ineoherent, hopping-type transport modeled by 
Levieh-Jortner rates, explieitly aeeounting for the inner and outer relaxation energies related to thermally 
indueed lattiee effeets and enhaneed eleetron-phonon interaetion at elevated temperatures. Complementary 
band-strueture ealeulations yielding temperature dependent effeetive mass tensors deviate stronger from ex¬ 
perimentally observed spatially anisotropie transport behavior. Thus, these results hint at the faet that by the 
partieular interplay of the transport energies the mobility of a given organie semieondueting material might 
appear to be band-like in a eertain temperature regime even though the underlying eharge earrier transport 
ean be of ineoherent, hopping-type nature. Building on this deseription, we further explore the role of the 
intermoleeular eleetronie eoupling and develop a proeedure to distinguish between its dependenee on statie 
and dynamie lattiee deformation upon temperature variation. 
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I. INTRODUCTION 


Charge carrier transport constitutes a decisive process in the operation of organic thin film 
devices such as light emitting diodes, photovoltaic cells and field-effect transistors (FETs).||T| De¬ 
spite tremendous efforts to develop high mobility organic semiconductor materials and their first 
commercial applications in displays and lightning, the microscopic mechanisms mediating the 
motion of electrons and holes on a molecular level and their interference with effects on spatially 
extendend length scales, in particular static or dynamic lattice deformations, have not yet been 
fully elaborated. A major challenge of up-to-date theoretical models is posed by the delicate en¬ 
ergetic balance between the individual transport processes, demanding for reliable mobility data 
over a broad temperature range and along different crystallographic directions to cope also with 
the spatial anisotropy of charge carrier motion in crystalline molecular solids. [O In this respect, 
organic single crystals constitute reference systems due to their material-inherent electronic behav¬ 
ior with only minor influences by chemical or structural inhomogeneities. [[3]| In combination with 
advanced simulations they enable detailed understanding of the charge carrier migration processes 
in relation to molecular packing, orientation and localization.® 0 A detailed overview dealing 
with the question of charge localization in different transport regimes has recently been provided 
by Fratini et. al.l0 

An elegant way to access these transport properties, that has recently been rendered suc¬ 
cessful, is based on the lamination of micrometer thick molecular crystals on vacuum-gap 
poly(dimethylsiloxane) (PDMS) substrates with predefined FET bottom-contact structures. 01 
To obtain the full information on the in-plane charge carrier mobility this experimental approach 
requires, in principle, just one crystalline sample without need for destructive readjustment.® 0 
But even though a band-like motion of electrons or holes might be anticipated and has been demon¬ 
strated for various polyaromatic single crystals e.g. by time-of-fiight measurements. ifTOl fTTl com¬ 
parable FET studies on laminated crystals predominantly reveal a thermally activated transport 
characterized by a mobility decrease towards lower temperatures. IfT^ This becomes even more 
remarkable considering the high in-plane mobilities of > 10 cm^/Vs observed already at room 
temperature in the case of laminated rubrene crystals. [IT3l The absence of a non-thermally acti¬ 
vated charge carrier motion renders it difficult to judge whether the observed behavior is caused 
by extrinsic obstructions such as morphologically or chemically induced trap states at the crystal 
surface or by intrinsic limitations, namely the small band dispersion and strong phonon coupling 
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expected for van-der-Waals bound solid states. [fT4ll To discriminate the various contributions to 
charge carrier migration and to advance existing transport models it is therefore mandatory to 
have organic semiconductor references showing a non-dispersive transport together with spatially 
anisotropic effects along the various crystallographic directions. 

We aimed for these challenges by investigating the temperature dependent charge carrier 
transport in A^W’-Bis-(heptafluorobutyl)-2,6-dichloro-l,4,5,8-naphthalene tetracarboxylic diimide 
(CI 2 -NDI) single crystals. This particular naphthalene diimide derivative has not only become a 
potential candidate for complementary logical circuits by its auspicious high electron mobility of 
1 cm^/Vs in thin films ifTSlfT^ and up to 8.6 cm^/Vs in solution processed micro-crystals. [fTTIl but 
can also be grown as single crystals by physical vapor deposition technique liTSl under ambient 
conditions [|T9ll as utilized in this paper. The habit of the sublimation grown CI 2 -NDI crystals is 
characterized by lateral dimensions along the (001) plane of several millimetres and thicknesses of 
several tens of microns making elastomeric FET measurements of the anisotropic in-plane mobil¬ 
ity feasible. Advantageously, both laminated crystal and PDMS substrate are expected to exhibit 
similar thermal expansion coefficients thus maintaining the structural integrity of the samples by 
avoiding significant interfacial strain at the examined temperatures between 175 K to 300 K. We 
restricted our transport studies to this range to avoid interferences with a polymorph occurring 
at lower temperatures. IIT^ Furthermore, operating the PDMS FET structures under vacuum as 
done in this study minimizes chemical surface reactions, which might lead to undesired effects by 
interfacial doping or dipole formation in the active channel region.® Altogether, these conditions 
allow for reliable investigation of the electron mobility which act as a test-bed for theoretical 
calculations of this key parameter as a function of temperature and crystallographic direction. 


II. EXPERIMENTAL 

Single crystals of the triclinic P-phase polymorph (PT, a = 5.346 A, b = 6.317 A, c = 19.157 
A, a = 92.63 °, p = 98.55 °,y= 109.00 ° at T = 300 K; CCDC: 1472770)|l20l were grown from 
purified CI 2 -NDI material by the physical vapor transport method [|T^ at 235 ± 5 °C under am¬ 
bient air. [fT9ll The transistor geometry utilized in this study consists of lithographically structured 
elastomeric PDMS supports covered by 3 nm thick Cr and 30 nm thick Au layers to form the 
contacts. By the elevated source and drain contacts a vacuum gap of about 4.9 ^um in height and 
channel lengths of 50 pm are generated. To improve charge carrier injection the pre-structured 
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FIG. 1. (a) CI 2 -NDI molecular structure and single crystal arrangement on top of an elastomeric FET 
support, (b) Molecular arrangement within the unit cell together with the relevant electron transfer directions 
a, b, a + b, a — b. 

Au contact pads have been surfaee funetionalized by spin eoating an aqueous solution of 0.004 
wt.% polyethylenimine ethoxylated (PEIE) onto the substrates. ll^ All /(E)-measurements have 
been earried-out in transistor geometry (i.e. with two probe souree and drain eontaets) under high 
vaeuum at a base pressure of 10^^ mbar in a Janis Researeh ST-500 mieromanipulated probe sta¬ 
tion using an Agilent 4155C parameter analyzer. Eleetron mobilities were ealeulated in the linear 
regime of the transfer eharaeteristies using Eq. 



( 1 ) 


The temperature dependent erystallographie strueture of the bulk samples was determined by 
a Bruker D8 Quest diffraetometer, whereas the eomplementary assignment of the erystallographie 
axes was aeeomplished by a Bruker D8 Diseover diffraetometer, both using eopper Ka radiation. 

Eigure [^a) pietures the moleeular strueture of CI 2 -NDI together with the alignment of a 
sublimation-grown crystal laminated on top of the pre-patterned PDMS support. The pie-shaped 
contaet strueture enables determination of the mobility along defined crystallographic directions 
at an angular resolution of 30 °. The main eharge transfer direetions within the two-dimensional 
briekwall-type plane of slip-staeked CI 2 -NDI molecules, i.e. a, b, a + b and a — b, are displayed 
in Eigure [^b). 
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III. THEORETICAL DECRIPTION BY LEVICH-JORTNER RATES 


To provide theoretical access to the experimental transport data, charge carrier motion in or¬ 
ganic semiconductors has to be divided into two limiting regimes depending on the degree of 
localization.If the electronic coupling J between the molecules is weak, charge carriers are local¬ 
ized and transport can be described by non-adiabatic hopping rates depending on J^. ll2^ In this 
case the time for the electronic wave function to move from donor to acceptor Xei = h/J (given 
by the reduced Planck constant h and the electronic coupling J) is much larger compared to the 
time for nuclear motion of the promoting mode Xnuc = 27t/(0. The high temperature limit of clas¬ 
sical charge-transfer theory assumes that all promoting modes are harmonic and can be treated 
classically yielding the often used Marcus rate, where the activation energy is mainly given by the 
reorganization energy ?i. ll^l2^ Here, the latter is a sum of two contributions: a larger one caused 
by fast relaxation of nuclear coordinates of the two involved molecular species upon changing 
their charge state (intramolecular and a smaller one originating from slow relaxation of the 
surrounding molecules (outer-sphere 

At ambient conditions, however, the fast intramolecular modes (typically C-C bond stretches 
with 0.2 eV 3> k^T) should be treated quantum-mechanically which is provided by the 

Levich-Jortner hopping rate. ll25l where the activation energy is determined mainly by 

2% |7|2 - S^exp(-5) ^ ^ | 

Here it is assumed that the donor is initially in its vibrational ground state, while N describes all 
accessible vibrational states of the acceptor at energy with S = being the Huang- 

Rhys factor, k^T the thermal energy and AE the site-energy difference. 

While the hopping picture is appropriate for disordered amorphous organic semiconductors 
with large energetic disorder and low average electronic coupling 7 1 meV, in organic crystals 

the strong n-n interactions can lead to 7 ~ 0.1 eV rendering it apriori not obvious whether a 

hopping model is applicable or not. Therefore, we have compared the highest electronic coupling 
of J = 82 meV, which occurs along the a-direction in our CI 2 -NDI crystal at room temperature, 
with the following localization energies: 

(i) intramolecular reorganization energy = 370 meV, describing the relaxation of nuclear 
coordinates of a molecule upon (dis)charging in vacuumllSll, 

(ii) outer-sphere reorganization energy = 40 meV, describing relaxation of the nuclei of 
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surrounding molecules in the crystal after the charge-transfer process [| 26 ] - I 2 ^ . 

(iii) difference in polarization energy AP = 300 meV for an electron localized on a single 
molecule within the molecular crystal versus spreading the charge across a dimer, 

(iv) energetic disorder a = 40 meV, due to molecular vibrations associated with a distribution 
of electrostatic and polarization energies, 

(v) intermolecular fluctuations (lattice phonons) leading to a broad distribution of J with a 
median of < 7 > = 80 meV and a standard deviation of 31 meV along the a-direction at room 
temperature (10% of molecular pairs have J < 50 meV and for 2% of the pairs J < 30 meV). 

Concerning point (v), it has been shown that the broadening of the J-distributions due to en¬ 
hanced intermolecular fluctuations at elevated temperatures might lead to stronger charge localiza¬ 
tion, which in case of rubrene single crystals can explain the decrease in mobility with increasing 
temperature. [0 [2^ However, the ratio between intramolecular reorganization energy and aver¬ 
age transfer integral at room temperature along the direction of highest mobility amounts to / 
<J> = 370 meV / 80 meV = 4.6 in our CI 2 -NDI crystals, whereas for rubrene / < 7 > = 159 
meV / 142 meV =1.1. This results in an a priori stronger localization of charge carriers in our 
system and thus favors its description by a hopping rate model. 

J and were calculated by density functional theory (DFT). All DFT calculations in this work 
were performed using the B3LYP functional and the def-TZVP basis set in the TURBOMOLE 
package [l30l except for J where the functional was changed to BP86 for performance reasons. 

The energy contributions (ii)-(iv) were obtained by the Thole model OTI based on atomistic 
charges. For the dimer configuration in (iii) we have equally placed 1/2 of the anionic atomic 
charges on each of the two molecules and then used the Thole model to evaluate the polarization 
response of the surrounding molecular crystal. 

Due to their fast time-scales, polarization interactions and intramolecular reorganization renor¬ 
malize the transfer integrals 7 to smaller values f. It has been shown that for pentacene polariza¬ 
tion leads to J'/J = 0.64. Il3^ With = 96 meV and a mean vibrational energy of h(p = OAleV 
for this compound, the intramolecular reorganization J'/J = exp[—X™^/{2h(p)] yields an addi¬ 
tional reduction by 0.75. li33l The intramolecular contribution has been already accounted for in 
the Levich-Jortner model by the Huang-Rhys factor (eq. [^. As in [[3^ the polarization effect 
has been estimated isotropic, this polarization effect has been neglected in our simulations of the 
anisotropic mobility tensor. 

If we assume a similar effect of polarization in our CI 2 -NDI samples but consider the larger 
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^int = 37 Q nieV at the same mean vibrational frequency, we arrive at a renormalized average 
coupling along the direction of highest mobility at room-temperature of J' = 80 meV * 0.64 * 
0.33 = 17 meV. By comparing this value to the energetic disorder of o = 40 meV the localization 
length of the charge carriers in units of intermolecular spacing L/a can be computed according to 
L/a = 8(o/7')^^ '^,[l33l which finally results in a localization length of 2.4 molecular units. Since 
already the non-normalized integrals represent broad distributions and along the b-direction the 
integrals are even smaller, the tendency to localize is even further enhanced. This leads us to the 
conclusion that in our CI 2 -NDI molecular crystals a localized picture and thus a hopping transport 
mechanism are valid. 

In order to complement this picture, we finally compare our temperature-dependent mobility- 
tensors based on the Levich-Jortner hopping rate to effective mass-tensors from a perdiodic DFT 
calculation of the band-structure (section [TV]). 


III.l. Molecular dynamics 

We have used molecular dynamics (MD) simulations in order to study how temperature affects 
the morphology of the CI 2 -NDI crystal and thereby the electronic coupling between the molecules 
in the (001) plane. A force field based on Optimized Potentials for Liquid Simulations (OPLS) 
parameters [l34ll was used with atomistic partial charges obtained from DFT while unknown in¬ 
tramolecular potential energy surfaces have been parametrized again using DFT OTI The only 
modification to the original OPLS parameters was a reduction of the van der Waals radius of the 
Cl atoms by 12%. When starting the MD simulation from the experimentally known crystal struc¬ 
ture at r = 200 K the equilibrated unit cell deviates by less than 1% in all six lattice parameters 
(a,Z7,c,a,P,Y) and we are also able to reproduce the experimentally observed thermal expansion 
of the crystal as shown in Figure As mentioned before, in this temperature range the crystal 
structure of sublimation grown CI 2 -NDI single crystals retain its triclinic symmetry [|T^ and thus 
differs from the monoclinic structure previously reported for solution grown crystals. ifTTII As indi¬ 
cated by Figurej^the experimental data for the thermal expansion coefficients, estimated by single 
crystal X-ray diffraction, can be well reproduced by our simulations based on molecular dynamics. 

MD simulations were performed using GROMACS [l35ll for boxes containing 18x18x6 CI 2 - 
NDI molecules with periodic boundary conditions. We simulated at constant particle number, 
pressure and temperature (NPT ensemble) with anisotropic Berendsen pressure coupling ll^ 
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FIG. 2. Relative lattice expansion AL/L{300K) along the crystallographic axes determined by X-ray analy¬ 
ses (full symbols) and force field simulations (open symbols) together with the respective thermal expansion 
coefficients, indicating the good agreement between theory and experiment within the examined tempera¬ 
ture range. 

(reference pressure 1 bar, time constant of X/> = 5 ps and compressibility 4.5 x 10^^ bar ^^) in 
combination with stochastic velocity rescaling for temperature coupling 071 (Xy = 0.07 ps) and a 
time step of 0.002 ps, using constrained bonds and the grid-based particle mesh Ewald technique 
for electrostatic interactions. 08l Morphologies for lower and higher temperatures were achieved 
by thermal annealing starting from the equilibrated MD structure at T = 200 K with maximum 
rates of 25 K/ns. By this protocol the deviation between equilibrated simulation and experiment 
for the most relevant lattice parameters within the transport plane a,b,yis less than 1% for r=175K 
and r=300K. 

III.2. Temperature dependent transport parameters 

In order to separate the different thermal effects on the anisotropic charge carrier transport 
we have, at first, utilized the lattice parameters from the MD-equilibrated crystal structure for 
r=300 K and r=175 K, keeping all molecules exactly at their lattice sites and with orientations 
defined by the crystallographic angles excluding any disorder ("frozen crystal"). To model the 
temperature and direction-dependent mobility we employ Eq. where we assume = 0.2 eV 
and compute all parameters for electron transport by ab — initio methods llSTIl yielding = 0.37 
eV (independent of temperature and direction) while all other transport parameters are summarized 
in Table m 

In the frozen crystal (labeled cry in Table |I]) the transport parameters are identical for all pairs 
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that belong to the same crystallographic direction shown in Figure [Jb). The small intermolecular 
distance d along the a-axis leads to the strongest electronic coupling J and lowest both fa¬ 
voring transport along this direction. Due to the nodal structure of the LUMO (lowest unoccupied 
molecular orbital) the coupling along a-\-bh negligible although d is not too large indicating that 
charge transport in this direction has to occur via sequential hopping in the a and Z?-direction. The 
coupling along the b- and a — Z7-direction is slightly weaker compared to that along a which effec¬ 
tively yields a two-dimensional transport within the brickwall-type packed CI 2 -NDI (001) plane. 
As the LUMO is located on the conjugated core of the molecule, transport along the c-direction, 
i.e. along the side chains, is neglected in our simulation. The biggest difference for the two prin¬ 
cipal transport directions a and b of the frozen crystal at 300 K and 175 K is the enlarged J along 
the a-direction due to the anisotropic contraction of the crystal upon cooling. 

As we have only one molecule per unit cell the energy difference vanishes, i.e. /S.E = 0. Fur¬ 
thermore and decreases with intermolecular distance J. OTII Note that 

as we have used a Pekar-factor of Cp = 0.05 OTl (which linearily scales such that only the 
first vibrational states contribute to Eq. 

The directional mobility iu={v-F)/F^ dX an angle (|) to the a-axis can now be obtained from all 
hopping rates performing in average 100 Monte-Carlo simulations for a single, randomly injected 
electron in a periodic box projecting its mean velocity v on the applied field ofF = 10^ V/m. 

For a one-dimensional transport along an arbitrary direction (z) Levich-Jortner-hopping along 
(k_|_) or against (k_) a small electric field relates to a small energy difference AF = eFAz —)■ 0. 
Implying that only N=0 significantly contributes in Eq. [^yields the following temperature depen¬ 
dence of mobility: 


b = 


{k+-k-)Az 

F 


27t:7^exp(—5) e{AzY f ^ 

hy/An'koutkBT k^T ^ \ ^k^T) 


niT) ~ 


exp(-£7) 

7 W 2 


( 3 ) 


Note that eq.j^also can be deduced from the Marcus rate model, however, with X = -I- ~ 

0.45 eV in the exponent a decrease in mobility by more than one order of magnitude will occur 
upon cooling, contradicting the experimental findings (see Figure Qa) and discussion further be¬ 
low). In the case of Levich-Jortner hopping, is small compared to AkgT (see Table |p. There¬ 
fore, the mobility is expected to follow the power law /r ~ 7^3/2 for a temperature independent 
fi. Note that this explanation of reduced mobility at higher temperature is different compared to 
the mechanism suggested in case of rubrene single crystals, where the increasing disorder due to 
stronger intermolecular fluctuations is assumed to decrease the localization length and thus the 
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direction 

r[K] 

d [A] 

J [meV] 


^E 

[meV] 


^out 

[meV] 



cry 

<•> 

a 

cry 

lg[72] o 

<lg[72]> 

a 

cry o 

<•> 

a 

cry 

<•> 

a 

a 

300 

5.35 

5.36 

1.71 

82 

-2.16 

-2.14 

0.37 

0 

0 

38 

41 

41 

2 


175 

5.20 

5.20 

1.12 

102 

-1.99 

-2.03 

0.21 

0 

0 

29 

38 

38 

2 

b 

300 

6.31 

6.30 

1.85 

58 

-2.47 

-2.43 

0.80 

0 

0 

35 

52 

51 

3 


175 

6.38 

6.38 

1.27 

57 

-2.48 

-2.18 

0.39 

0 

0 

25 

53 

51 

2 

a — b 

300 

9.45 

9.45 

1.53 

16 

-3.55 

-3.70 

0.43 

0 

0 

40 

91 

90 

2 


175 

9.40 

9.41 

1.03 

18 

-3.49 

-3.56 

0.24 

0 

0 

30 

79 

78 

2 

a + b 

300 

6.90 

6.90 

1.45 

0.2 

-7.51 

-7.16 

0.96 

0 

0 

39 

81 

80 

2 


175 

6.87 

6.85 

1.03 

0.1 

-8.05 

-7.30 

0.48 

0 

0 

29 

88 

89 

2 


TABLE I. Transport parameters for Eq.|^in the frozen crystal (cry) as compared to median values <•> and 
standard deviations a for the distributions induced by lattice phonons along crystal directions defined in 
Eigure [TJb) for different temperatures. Distributions of center of mass distance d, site-energy difference 
AE and outer sphere reorganization energy E™* are Gaussian. As electronic couplings J are very broad and 
non-Gaussian we display the median and a for logjQ[(J/eV)^]. Thus less negative numbers correspond to 
higher J. Eor comparison we have also converted J for the frozen crystal into logjo[(J/eV)^] with largest 
deviations from the frozen crystal occurring at 175 K along the b-direction. 

charge carrier mobility. O But the predicted T ^ dependence of the mobility is not supported by 
the experimentally observed behavior of our samples. In contrast, the Levich-Jortner rate 
model can reproduce the measured /j(r) data and therefore, appears to be better suited. 

Finally, we also account for lattice vibrations (phonons): Our MD-simulation reveals that all 
molecules undergo vibrations around their lattice points which are well described by Gaussian 
distributions of center-of mass distances d (median and standard deviation are labeled <•> and a 
in Table |I} respectively). 

As a consequence of these phonons each molecular pair acquires a unique conformation and en¬ 
vironment which transform the anisotropic transport parameters into distributions. We have taken 
two representative snapshots, each containing 6x18x18 molecules, from our MD simulations after 
equilibration at 175 K and 300 K, respectively. Then all Jortner rates for our KMC simulations 
were evaluated by the parameters obtained from the molecular pairs within these snapshots. Site 
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energy differenees AE eomputed ineluding eleetrostatie and polarization interaetions [|3T]| are well 
deseribed by Gaussian distributions (Table with Oa£ representing the energetie disorder. As the 
energetie disorder is similar for all direetions, its inelusion in Eq. [^does not lead to a rotation or 
deformation of the tensor. Moleeular vibrations also result in Gaussian broadening of due to 
different environments but their effeet is weak as the standard deviation is small and isotropie. 

While the dependenee of eleetrostatie and polarization effeets (both affeeting AE and on 
intermoleeular distanees is given by a power law, the eleetronie eouplings J depend exponentially 
on distanee. Furthermore, via the nodal strueture of the orbitals J is also very sensitive to the 
relative orientation between two moleeules. Therefore, the main influenee on shape and anisotropy 
of the mobility tensor results from the temperature dependenee of the J distributions. In ease of 
CI 2 -NDI, just ineluding distributions in J on top of the statie pieture (while keeping AE = 0 and 
fixed) does not lead to strong ehanges in absolute mobilities whieh is indieated by the small 
aspeet ratio of ^ around 3 (Figure |^) allowing for effeetive 2 dimensional transport within the 
(001) plane. As a eonsequenee, the few weak eouplings in the low energy tail of /-distributions 
do not result in bottle-neeks for eharge transport. 


IV. BAND-STRUCTURE CALCULATIONS 

Eleetronie-band-strueture ealeulations for the two erystallographie struetures obtained exper¬ 
imentally at 200 K and 300 K were earried out using DFT and the Perdew Burke Emzerhof 
Generalized Gradient Approximation GGA-PBE [|^ l40l as implemented in Atomistix Tool Kit 
tATKl. [l4Tl The k-point sampling was 6x6x2 for both struetures, and we set 75 Hartree as the eut- 
off energy and seleet 10^^ Hartree as the eonvergenee eriterion for the total energy. The eleetron 
temperature is taken as 300 K. Upon obtaining the eleetronie-band-strueture for both struetures, 
the fraetional k-point X(0.5, 0.0, 0.0) was identified as a eonduetion band minimum. For the eal- 
eulation of the effeetive mass tensor and the effeetive mass direetion, ATK’s effeetive mass tool as 
implemented in version 2014.2 was used. 


V, EXPERIMENTAL MOBILITY TENSOR 

Performing the eleetrieal eharaeterization in FET geometry and in the temperature range de¬ 
seribed above. Figures [^a,b) show representative transfer and output eharaeteristies of a CI 2 -NDI 
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FIG. 3. (a) Transfer characteristics measured at 175 K and 300 K along 253 ° (low mobility) and 343 ° (high 
mobility) with respect to the a-axis at Fd = 1-5 V, and (b) corresponding output characteristics measured at 
Fg = 5 V. 


single crystal FET recorded at 175 K and 300 K along the directions of low (253 °) and high 
(343 °) mobility (angles are always declared with respect to the a-axis). As indicated by the small 
hystereses in the /(F)-curves upon up- and down-cycling, dielectric relaxation losses caused e.g. 
by trap-assisted charge carrier localization or pronounced injection barriers are of negligible rel¬ 
evance for the deduced transport parameters. For single transistors, electron mobilities estimated 
in the linear transfer regime at 175 K by means of Eq. [grange between 0.3 cm^/Vs to 2.8 cm^/Vs 
depending on the respective crystallographic direction. The well-defined On/Off states in our sin¬ 
gle crystal FETs are confirmed by their related current ratio of about 10^ which can be attributed 
to the small leakage currents of the utilized vacuum gap structure. Threshold voltages of about 3.5 
± 0.6 V indicate the presence of trap states which can be effectively filled already at a low density 
of accumulated charges considering the small gate-voltage and low capacitance of 0.18 nF/cm^. 
Accordingly, at 300 K mean threshold voltages are reduced to 2.1 ± 0.6 V mainly due to thermal 
detrapping of electrons. The activation of otherwise immobili z ed charges at elevated temperatures 
is confirmed by enhanced leakage currents at room temperature leading to On/Off ratios of about 
ten and by a superimposed ohmic contribution in the output characteristics. Most strikingly, at 
300 K the lowest and the highest in-plane mobility are reduced to 0.1 cm^/Vs and 1.5 cm^A'^s, 
respectively. 

Measuring the FET characteristics along various directions in the CI 2 -NDI (001) plane the 
electron mobility tensor has been constructed. Exemplarily, the anisotropic electron mobility at 
300 K and 175 K is displayed in Figure Qa). To estimate the range of confidence, mobility values 
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FIG. 4. (a) Electron mobility tensors of CI 2 -NDI single crystals measured at 300 K (magenta) and 175 K 
(cyan), (b) The corresponding simulated tensors are calculated by a hopping-like transport model based 
on Levich-Jortner rates. The simulation includes static as well as dynamic lattice effects (see below: sec¬ 
tion VII|). (c) Effective mass-tensors from periodic DET band-structure calculations based on experimental 


crystal structures at 300 K and 200 K. 


along different direetions have been averaged in eomplianee with the PI point group symmetry 
of the CI 2 -NDI P-phase erystal structure (see above), with deviations visualized by the error bars. 
At room temperature the angular dependent mobility resembles a constricted ellipse with its long 
axis tilted by -22 ° and its short axis by -112 ° with respect to the a-axis. The mobility aspect 
ratio along the principal directions A^max//rmin amounts to three and thus resembles that of other 
organic compounds like Pentacene. lH^ Rubrenell9l[T3l or DCNQ-Q Upon cooling to 175 K, the 
overall electron mobility increases by a factor of two with the principal axes of the mobility tensor 
remaining fixed within the angular resolution of ± 5 °. 


VI. SIMULATED MOBILITY TENSOR 

Complementary to the experimental results, the anisotropic mobility tensor at room temperature 
has been calculated by the theoretical Levich-Jortner rate model introduced in detail in section [IlH 
The predicted electron mobilities as well as their spatial anisotropy are illustrated in Figure |^b) 
exhibiting a good agreement with the experimental data, in particular concerning the difference in 
the absolute mobility values by only a factor of two. At room temperature both experiment and 
simulation reveal the direction of highest mobility at an angle of 338 ° and display a pronounced 
constriction in the direction of lowest mobility. At 175 K the overall mobility values increase by a 
factor of two leaving the aspect ratio Hm&x/Pmm almost unchanged, except for a partial lifting of the 
constriction along the direction of lowest mobility. For all of these transport characteristics, which 
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300 K 

175 K 

175 K 

200 K 

Change 


Exp 

Hopping 

Band 

Exp 

Hopping 

Band 

Exp 

Hopping 

Band 

A^max (cm^A/s) 

1.2 

2.8 


2.4 

4.8 


x2.0 

X 1.7 


‘l)max (°) 

158/338 

158/338 

161/341 

162/342 

152/332 

158/338 

+ 4 

-6 

-3 

Amin (Cm2/Vs) 

0.4 

1.1 


0.7 

2.3 


X 1.8 

X 2.1 


(°) 

68 / 248 

68 / 248 

71/251 

72 / 252 

62 / 242 

68 / 248 

+ 4 

-6 

-3 

A^max / A^min 

3.0 

2.5 

1.6 

3.4 

2.1 

1.7 

X 1.1 

xO.8 

X 1.1 


TABLE II. Measured and calculated mean values of the lowest and highest mobility as well as of their 
angular relation with respect to the crystallographic a-direction at 300 K and 175 K (band: 200 K). The 
variations upon cooling are listed in right column and indicate the good agreement between experiment 
and theoretical predictions by the Levich-Jortner rate model except for some slight deviation in the tensor 
rotation at low temperatures. 

are summarized in Table the good agreement between the experimental data and the simulations 
eorroborate the suitability of the Levieh-Jortner rate model. Note that the use of Mareus rates 
would lead to deereasing mobilities with deereasing temperature without affeeting the rotation of 
the tensor signifieantly. Therefore, the temperature dependent behavior of the eleetron mobility 
tensor establishes fundamental eriteria to judge on suited transport models and thereby, on the 
essential mieroseopie eharge migration meehanisms. Within the ehosen approaeh we were able to 
prediet these variations, whieh are mainly determined by J^. 

In order to eomplete the pieture we have performed band-strueture ealeulations by periodie 
DFT based on the experimental erystal struetures at 300 K and 200 K (details are given in See- 
tion sg. From the band-strueture we extraeted the inverse effeetive eleetron mass-tensor, whieh 
is displayed in Figure |^e). In band-transport the mobility ^ = ex/m* depends on the eleetron 
eharge e, the seattering time x and the inverse effeetive mass 1/m*. Beeause x ean eover several 
orders of magnitude [l43ll we deeided to not ealeulate the mobility. Instead, one ean eompare the 
direetion and shape of the inverse effeetive mass tensor to the mobility tensors in Figure 
if one assumes x to be isotropie and temperature independent. At 300 K the direetion of highest 
mobility is found at an angle of 341 °, the direetion of lowest mobility at 251 °. In eomparison to 
the simulations based on Levieh-Jortner rates only a small differenee of 3 ° is observed for both 
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directions, which shows that the electronic coupling is related to the band structure. This match 
is to be expected as both models rely on fi. Going towards 200 K the directions of highest and 
lowest mobility as well as the aspect ratio do not change significantly. 

The biggest difference between the inverse mass tensor and the mobility tensor from the Jomter- 
rates is (i) a less-pronounced constriction along the a 4- h-direction especially at 300K and (ii) a 
lower aspect ratio i^max/l^min in case of the former (see Table |^. As the mobility tensor deduced 
from the Jortner-rates qualitatively agrees better with the experimental findings in those two points 
we assume a localization picture to be justified. 

However, despite the remarkable agreement between experimental and Levich-Jortner mobility 
values at 300 K and 175 K and their overall spatial anisotropy, a small but nevertheless important 
deviation can be realized by the rotation of the principal axes upon cooling which shed light on the 
underlying transport processes. While in the measurements the tensorial rotation at 175 K occurs 
to be small (< 4 ° towards the a-axis) and within the experimental error, the theoretical calculations 
yielded a rotation by about 6 ° and, even more strikingly, away from the a-axis. As this effect 
appears to be of thermal origin, we will address in the following the two major contributions to 
the temperature dependent transport in organic single crystals, namely thermally-induced lattice 
distortions leading to a static variation in the LUMO overlap along the anisotropic intermolecular 
distances within the unit cell [[44]| and, secondly, the thermal occupation of phonon modes and 
their spatially anisotropic coupling to charge carrier motion. [|4l 


VII. STATIC AND DYNAMIC PICTURE 

To judge on the first contribution to the temperature dependent mobility, we investigated the 
thermally induced lattice expansion, i.e. the static picture, along the main crystallographic direc¬ 
tions of the CI 2 -NDI unit cell between 100 K and 300 K as described in section [IIlT] For this 
purpose, the well-known crystal structure at T = 200 K was taken as a starting point for further 
equilibration using a force field based on OPLS parameters. 1(341 The converging parameters of the 
unit cell coincide with the experimental ones within an error of 1 %. From here, all other simulated 
crystal structures were obtained by thermal annealing. Disregarding fluctuations in inter-molecular 
distances and orientations (phonons) in the static picture (static LUMO overlap), we could inves¬ 
tigate the structural influence on the charge transport separately. As key result, upon cooling from 
300 K to 175 K the mobility tensor rotates by 8 ° towards the a-axis (see Figure]^ a)), which can 
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be explained by the contraction of the a-axis and the simultaneous expansion of the h-axis upon 
cooling (see Figure]^, which results in an enhanced electronic coupling in the a-drrection com¬ 
pared to the Z7-direction at 175 K, as can be seen in Figurej^b). For this purpose, we have indicated 
the coupling elements of the static picture by lines in the lower part of this figure as compared to 
their distributions in presence of phonon modes above (which will be discussed below). We have 
omitted the a -\- Z7-direction as by the nodal structure of the LUMO in combination with the face- 
to-face orientation of molecular pairs the squared coupling is on average five orders of magnitude 
weaker than e.g. for the a-axis. Due to the contraction along the a-axis at low temperatures, the 
static picture results in improved coupling along the a-direction yielding a rotation of the mobility 
tensor towards this axis. 

Mobilities are larger in the static case (Figure |^a)) compared to the dynamic case of (Fig¬ 
ure Qb)) due to energetic disorder which is missing in the former (see section III.2| ). The increase 
of the mobility upon cooling by a factor of 2.4 to ^Umax = 17.6 cm^/Vs can only partly be explained 
by the 1.5-fold increase of the rate-relevant 7^ along the a-axis. The missing contribution can be 
deduced from the temperature dependence of Eq. This Equation indicates that changing T from 
300 K to 175 K will increase /r along the a-direction by /J = 1.7 due to small = 0.04 eV, 
while along the a — h-direction the increase by = 1.1 is weaker due to larger = 0.1 eV. 
Combining the gain in along the a-direction (//^=1.5) (see Table|^ with the thermally induced 
gain in the rate (/J = 1.7) we are able to explain the overall increase of /Umax by 2.4. This sig¬ 
nificant increase along the a-direction compared to the h-direction and the resulting rotation of 
the mobility tensor in the static picture is not reflected in the experimental tensor. Therefore, the 
static lattice deformation alone is not sufficient for the description of the electronic coupling upon 
temperature variation. 

Now we adress the second effect, the presence of phonon modes, i.e. thermally induced fluc¬ 
tuations of inter-molecular distances and orientations. They broaden the delta-like shape of the 
electronic coupling elements which have sufficed the static picture accounting solely for thermal 
expansion. Therefore, the activated phonon modes lead to a reduction of mobility (see Eigurej^b)) 
by a factor of 2.6 for T = 300 K and by 3.7 for T = 175 K with respect to the mobility of the 
frozen crystal approach (Eigure|^a)). The reduction is stronger in case of lower temperature as 
6 = (^aeH^bT) = 2.0 and 1.6 for T = 175 K and 300 K, respectively, and jj ~ exp[—(Cd)^] [|45ll 
with a constant C of the order of one. 

The distributions of J are non-Gaussian and very broad as shown logarithmically for the rate- 
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FIG. 5. (a) Simulated mobility tensors using morphologies from experimental crystal structures at differ¬ 
ent temperature, however without molecular vibrations, (b) Distribution of electronic coupling elements 
in a static picture (lower part) and in presence of non-local phonon modes (upper part) together with their 
temperature dependent redistribution along different directions. While at 175 K the static picture shows 
improved coupling mainly along a, including phonon modes also improves coupling along b due to a redis¬ 
tribution of weight from the low-energy tail as indicated by the median shift of the distributions displayed 
by the upper lines and arrows. 

relevant in Figure [^b). For the prineipal transport direetions a, b and a — b, the biggest effeet 
upon cooling occurs along b, where the distribution of J is broadest: the median is shifted to much 
higher values and the low-energy tail is reduced significantly upon cooling although the mean of 
increases only slightly (and therefore Jh in the frozen crystal does not change upon cooling). The 
final mobility tensors including all vibrational broadenings related to J, AE, and are displayed 
in Figure Qb) and show that the rotation towards the a-axis caused by thermal contraction of the 
crystal at lower temperatures (Figure |^a)) is outbalanced by the changes in distributions of J 
which leads to an overall rotation away from the a-axis. 

Comparing the theoretical predictions on the temperature induced static effects alone (expan¬ 
sion of the lattice —)■ rotation by 8 ° towards a-axis) versus those including also dynamic effects 
(expansion and phonon modes —rotation by 6 ° away from the a-axis) one rationalizes that the 
two effects compete against each other with the resulting rotation of the tensorial mobility at low 
temperatures depending on their respective size and thus being characteristic for the organic ma¬ 
terial under study. In particular, for the anisotropic transport behavior observed in our CI 2 -NDI 
single crystals we identified almost a compensation of both effects leading to a spatially fixed 
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alignment of the prineipal axes of the eleetron mobility tensor (Figure [Ja)), whereas the simula¬ 
tions hint at a slight overestimation of the broadening of the eleetronie eoupling distribution by 
phonons and its eontribution to the eleetronie transport (Figure |^b)). 


VIII. SPATIAL ANISOTROPY OF TEMPERATURE DEPENDENT MOBILITY 


Finally we adress the temperature dependent mobility along the prineipal erystallographie axes 
by fitting with the Hoesterey-Letson (HL) model. [|4^ This model implies the thermal release 
of eharge earriers being loealized by Boltzmann-distributed traps at an effeetive energy and 
density A^r- Even though the HL-model was introdueed originally for oligoaeene single erystals 
showing a band-like transport, the intrinsie mobility in this deseription, can be adapted to the 
respective conduction behavior, i.e. to either band- or hopping-type transport, with its temperature 
dependence described by a power law. This law typically accounts for the dominant scattering 
mechanism in band transport by its exponent n. WH In this case, we expect 1.5 in eq. Q not 
because of the scattering by acoustic phonons but for the temperature dependence of the Levich- 
Jortner model in eq. [3]for the case that < Ak^T. The denominator limits the mobility at low 

temperatures due to increased trapping times 




1^0 ■ 


^ + Nr 


(4) 


Figure [^a) displays the measured temperature dependent electron mobility, which can be rea¬ 
sonably fitted by the described model and reveals only the transport along the direction 43 °, i.e. 
close to the a -I- A-direction, to be limited at lower temperatures. The latter is caused by either trap 
states or the poor electronic coupling along this direction (see above). The obtained temperature 
exponent n shows a distinct relation to the crystallographic transport directions as indicated by its 
angular variation in Figure [^b). The value of n ~ 1.5 along the a-direction and n 1.3 along 
the h-drrection assorts well with the expected value of n ~ 1.5 according to eq. The reduced 
value in the A-direction might stem from the increased which weakens the temperature depen¬ 
dency of the mobility due to the exponential term in eq. Along the a -|- h-direction the exponent 
n is about twice as large. As the electronic coupling J is orders of magnitudes lower between 
molecules along this direction (see Table charge carriers preferentially perform two consecutive 
hops along the a- and h-direction, which might lead to the increased value of n. 
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FIG. 6. (a) Temperature dependent electron mobility measured along 163 103 73 ° and 43 ° with respect 

to the a-axis and fitted by the Hoesterey-Letson transport model. Below 200 K trap limited transport occurs 
only at 43 i.e. close to the a + b direction, (b) Absolute values of the power law exponent n along different 
transport directions. To illustrate the entire angular behavior, half of the data points is mirrored according 
to the inversion symmetry of the CI 2 -NDI single crystals. 

IX. CONCLUSIONS 

We presented a eombination of approaehes on the temperature-dependent anisotropic electron 
mobility in CI 2 -NDI single-crystals by means of FET measurements and theoretical modeling 
of the microscopic transport via the Levich-Jortner rate model. These simulations demonstrated 
that the rotation of the mobility tensor with respect to the crystallographic orientation provides a 
criterion on the thermally induced contributions by static and dynamic lattice deformations. We 
have also compared the experimental results to band-transport from periodic DFT calculations 
and find a better agreement for the hopping approach. As a main result, we were able to show that 
hopping can be the dominating charge transfer mechanism despite the observed dependence 
of the mobility. Thus, our results indicate that for a given organic semiconducting material an 
increase of mobility at lower temperatures can also be explaind by an incoherent motion of charge 
carriers thoroughly taking into account static and dynamic lattice effects. 
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